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Abstract 

We point out that, due to the use of static nucleon propagators in Heavy Baryon Chiral Pertur- 
bation Theory, the current calculations of the nucleon-nucleon potential miss certain contributions 
starting at two loops. These contributions give rise to contact interactions, which are both paramet- 
rically and numerically more important than the so called NNLO potentials. They show a peculiar 
dependence on the light quark masses, which should be taken into account when performing chiral 
extrapolations of lattice data. However, they do not appear to have an impact on phenomenology 
since they can be absorbed into redefinitions of unknown parameters which are usually fitted to 
data. 
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I. INTRODUCTION 



Since Weinberg's pioneering work in the early 90's 111, [M, there has been an enormous 
development of effective theory methods for few nucleon systems (see |3j, |4J, [5j for recent 
reviews). However, there is still a lack of consensus on how calculations in t 
effective theory including pio ns must be organized Jf], E, |[ 
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28|,l29|,l30j|3l|, 




391 ] . Starting 

from the Heavy Baryon Chiral Lagrangian (HB^L) 40|, one is interested in building a low 
energy effective theory for nucleon energies much lower than the pion mass. It seems clear 
that the remaining effective theory must consist of nucleons interacting through a potential. 
Hence the program may be divided in two: (i) calculating the potential and (ii) organizing 
the remaining quantum mechanical calculation. We shall not enter in point (ii) here, which 
is where the main difficulties arise. Concerning point (i), following the so called Weinberg 
approach, it is commonly believed that one can use HB;\PT counting rules and the outcome 
of the calculation can be easily organized in standard chiral counting in 1/A X ; A x ~ 47iF n , 
F n being the pion decay constant. We point out here that this is not so. If we understand the 
potentials as matching coefficients which arise after integrating out higher energy degrees of 



freedom 
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411 ] . then, starting at two loops, there are contributions to them which are missed 



if the static approximation is used for the nucleon propagators, as prescribed by HB^PT 
counting rules 62] . This is due to the fact that the energy scale given by the pion mass 



has an associated three-momentum scale ^m^M^-, being the nucleon mass, which may 
make the kinetic term of the nucleon prop agator as important as the energy when integrating 
out degrees of freedom of energy 63J]. This becomes apparent when individual Feynman 
dia gram s are analyzed with the method of the strategy of regions or threshold expansions 



45|, 



46] with non-relativistic nucleon propagators, as will be shown in the next section. This 
kind of contributions are related to the so called radiation pions discussed in Refs. [15j, \l& 

nn 

within the KSW approach|13l. Il4j. They are also related to the fractional powers of the pion 
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521. 



mass which appear in pion production off nucleon-nucleon systems 
However, to our knowledge, they have never been discussed in the framework of the original 

. We present in section HTT1 their complete 
contribution at two loops. In order to obtain it, the leading order lagrangian (in the isopin 
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limit) augmented by the kinetic terms of the nucleon suffices, 



(1) 



The purely pionic sector reads 

C w = S- {TrfV^VU + n£(U + Ut)]} , U = e 4 ^ , (2) 

7r a is the pion field and r a the (isospin) Pauli matrices. The pion-nucleon sector reads 

£, w = JV*(iD„-^(u.f) + ^-)jV. (3) 

N is the nucleon field, gA is the axial vector coupling constant of the nucleon, <r the (spin) 
Pauli matrices, and, defining u 2 (x) = U(x), the covariant derivative of the nucleon field is 
given by 

D,N= (dp + ^^uljN , (4) 
and the axial- vector type object u M by, 

= i(u'Vfj,u — uV^) = i{u>, V^w} = iu^V fJJw . (5) 

Finally, the nucleon-nucleon sector reads 

C NN = -—N^NN^N - —rttrNrttrN (6) 

where Cs and CV are low energy constants to be determined from the experiment or from 



QCD, for instance by lattice simulations 53j, |54j]. In section HVl we discuss the relevance of 



these contributions in higher order calculations and the significance of our results. Section 
IVl is devoted to the conclusions. Details of our calculation are presented in the Appendix. 

II. THRESHOLD EXPANSIONS: A SAMPLE CALCULATION 

In this section we will introduce the method of calculation that we will apply. Many 
processes in physics, like those involving heavy quarks or our pion-nucleon interaction, involve 
more than one mass scale. Such processes are notoriously difficult to calculate in perturbation 
theory beyond the one- loop level. To proceed one has to resort to approximations, either 



numerical or analytical. Among the latter finds its place the strategy of regions 45|, |46(. The 



idea is to perform an asymptotic expansion of the integrals in certain ratios of mass scales, 
so that the resulting integrals appearing in the calculation of every term in the expansion 
are simpler, and the expansion is homogeneous (each integral appearing in the construction 
contributes to a determinate power in the expansion parameter). In short, the method goes 
as follows: 

1. Determine the large and small scales in the problem. 

2. Introduce factorization scales and divide the loop integration domain into regions 
in which the loop momentum is of the order of one of the scales in the problem. 

3. Perform, in every region, a Taylor expansion in the parameters which are small in 
the given region, and stay at leading order. At this point, keep only the relevant 
regions and discard the rest: the relevant regions are those that somehow maintain the 
structure of poles of the original integral. If we integrate over multiple momenta, and 
the integrand has several propagators, so that we have one or more poles associated 
with each momentum, after performing the Taylor expansion we should still have at 
least one pole for each momentum. If we end up loosing all the poles that were 
associated with one of the variables of integration, the region we are considering is 
irrelevant and must not be taken into account. 

4. After expansion, ignore all factorization scales and integrate over the entire loop inte- 
gration domain in every relevant region. 

The non-trivial point to justify is 4, which also guarantees the homogeneity of the expansion 
formula. In order for that point to hold it is essential to use dimensional (or analytic) 
regularization for the integral, even if it is finite in four dimensions. Loosely speaking, 
4 follows in dimensional regularization from the property that all integrals without scale 
vanish, but the truth is that at present day there are no mathematical proofs of the method 
of regions. The best we can say is that it has not failed yet, giving asymptotic expansions for 
any diagram in any limit and having been checked in numerous examples when comparing 
results of expansion with existing explicit analytical results. 
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We will now show how this method works and how it enables us to find new contributions 
to the nucleon-nucleon potential. Consider the following diagram, 




1 1 1 



(k — I + p — p 1 ) 2 — m 2 + ie I 2 — m 2 + ie k 2 — m 2 + ie fco _|_ „o _ ( k +p) 2 



P 2M N + lt 



where we have put ourselves in the center of mass frame, that is, 



ie 



(7) 



Pl = (P°,P) P3=(P°,fi) , c . 

5 5 {°J 

P2 = (p°, -p) Pa = (p°, -P f ) 

and so we write Pi ~ p-s = (0,p — f?) as just p — fi- We will use this notation whenever the 
time component of a four-vector is zero or neglected. 

We will assume that the nucleon momentum p and p', and the momentum transfer p — p' 
are of order whereas the nucleon energy p is of order m 2 /M^, which fixes the scales 
of the effective theory. The relevant energy scales in the diagram are m 2 /M^ and m n and 



the associated three-momentum scales and y/m^MN respectively. We will have in mind 
the philosophy described in {39]: modes with energies larger than m 2 /M N are integrated out 
giving rise to the potential. The strategy of regions will help us to separate these modes 



from the ones which must be kept in the effective theory (see 55j for an example of such a 



calculation). Now we must break the integration domain into pieces. The relevant regions are 
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easily found by inspecting the poles of the propagators: they are the parts of the integration 
domain in which one or more of the propagators inside the integral develop a pole. Let us 
display next the leading term of the Taylor expansion of the integrand in several regions: 



A; ~ 777^ , k ~ m„ 



r ~ , I ~ m„ 



[k — I + p — p'Y(k — I + p — p') 
{k — I + p — p 1 ) 2 — ml + ie 
1 1 



k° - 1° ~ m, 



fc'fc* 



P — m\ + ie k 2 — m 2 + ie 
1 1 



(9) 



(10) 



x- 



k° + ie k°-l° + ie k° — iek° — l° — ie ' 
This is a pure three pion exchange potential contribution, the same we find in the static 
approximation. Strictly speaking, one should make sense of the pinch singularities 
above. There are several ways to treat them in the literature (see for instance £ 



56 



57|]). In the EFT framework we are working in they have a natural interpretation, 



since the same singularities appear in the lower ener gy EFT, which arise from the 



56 



5811 . Once these iterations 



iterations of lower order potentials in the static limit 
are subtracted, a well defined expression is obtained. In this particular case, the 
needed subtractions are the iterations of an OPE (one pion exchange)- TPE (two pion 
exchange), TPE-OPE and OPE-OPE-OPE. The same comment holds for the various 
pinch singularities we will encounter below. 



k° ~ m w , k ~ m, 



1° ~ m n , / ~ m n 



(k-l + p- p')\k -l + p- pj 
— (k — I + p — p 1 ) 2 — m 2 
1 1 



- f ^ JUL 

h 1 M N 



IW 



k l k l 



I 2 — m 2 + ie k 2 — m 2 + ie 



1 



k* + ie k°-ie k o_io +p o_(j2zj+z£ + ie fc o _ p _ p o + <k£&? 



(11) 



(12) 



16 



2M N " r 1 2M N 

We recognize here the contribution of a two pion exchange potential followed by a one 
pion exchange potential, a contribution that also exists in the effective theory (i.e. in 



6 



the quantum mechanical calculation) and hence, to avoid double counting, must be 
discarded. 



1° ~ 171^ , I ~ 

(k-l+p-p'Yik-l+p-p'Y VI s 



k l k l 



X ■ 



(k-l+p-f) 2 -ml + ie l 2 -ml + ie -k 2 - m 2 
1 111 



(13) 



(14) 



k° + P° ~ i W 1 + ^ k°-p° + 



— ie 



1° 



ie 



ie 



2M N r 1 2M N 

As opposed to the former region, this one gives a one pion potential followed by a two 
pion potential, and must be discarded likewise. 



(15) 



70 r 



(k - l + p - p')\k -l + p- p')" 
— (k — I + p — pi) 2 — m 2 

1 



VI s 



lut 



k l k 



_/2 _ m 2 _ k 2 _ m 2 k +p 0_ (|+gi + le 



X 



2M N 



2M N 



(16) 



This is a one pion potential iterated three times, which also exists in the effective 
theory, and hence is to be dropped like the other contributions. 



k° ~ m w , k ~ ^/M^vtlt, 



k° — 1° ~ m n , k — I ~ \[M^m % 



(k-lY(k-l) r VI s k l k l 1 

-(£-f) 2 -f 2 -fc 2 fc0__|_ + ; e A;0 + 



fc 2 



;i7) 



(18) 



7 



X- 



UO _ 70 _ i 7 V i-0 _ 70 i (^-0 2 _ ,V 

After performing the integrations over A; and 1° we are left with 

(k-iy(k-iy in* k l k i . _„ 1 



M 



A'" 



(19) 



-(A; - 1) 2 -/ 2 -A; 2 "k 2 (k-l) 2 

The remaining integrand has no scales, and therefore vanishes when we integrate over 
the remaining d — 1 dimensions for k and / using dimensional regularization. 

By analogous arguments, we find that all the remaining regions have also vanishing contri- 
butions, except for the following one. 



k° ~ , k ~ \fM N m lt 



k°-l° ~ 777, 



A;' A;* 



-A; 2 / 2 -m 2 +ie _p fc o _ fc^ + ie fc o + 



2M N 



k 2 
2M N 



(20) 



(21) 



— le 



x ■ 



k° - 1° - 



k 2 



+ ie k° -1° + 



k 2 



— le 



2M N 1 2M N 

This contribution, that is not the composition of one or two pion potentials, nor does it 
vanish, is the new piece of the nucleon-nucleon potential arising from considering non- 
relativistic (rather than static) nucleon propagators. Note that the key feature of this 
contribution is that one virtual energy (Z°) gets the scale from one of the loops and 

— * 

feeds it into the second loop so that the three-momentum of the latter (A;) gets in turn 



the scale \/M^m v . This large momentum scale makes the pion exchange interaction in 
that loop instantaneous. This is important in order to identify the relevant diagrams 
displayed in the Appendix, and also explains why at least two loops are necessary to 
see this kind of contributions. 

So, summarizing, after the analysis of poles of the propagator we have found a purely 
three pion exchange potential (which is the result of the static case), various iterations of 
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one and two pion exchange potentials, and a new contribution given by 

(2ir) d I 2 - ml + it 

(22) 



\F n J y 1 2 2 2 2 2 2 7 (27r) d 7 



1 1 



The tensorial part of our integral reduces to 



a\ o\ a\ a r 2 a{ a r 2 


a\ a[ o\ a\ a J 2 a l 2 a\ a{ a[ a l 2 o\ o\\ 


1 k 4 ! 2 


2 2 2 2 2 2 


222222 222222) 


1 (d-l) 2 (d + l) 



1 (23) 



d-1 2 6 



And so (122]) becomes 



F n J v 64 J (27r) d - 1 7 (27r) d ~ 1 

Mjk 1 1 

x —j= , • (24) 

fc8 ^f 2 + m2 fc2 + M N \l 2 + m 2 



Integrating over the solid angle 

/ = 

J.U//2) 

and defining 



2vr d / 2 

dn d = tttttt^: , (25) 



- M ^ P ! m ' and ,= * (26) 



we can transform the integrals into the product of two beta functions, and so we find that 
the new contribution is 



d-1 3d-7 

D 



(7ri ■ r 2 - 6) 



2 



r(¥) 

Regularizing in d = 4 — 2e dimensions this is 



Fj y 'd-l 64 2 2d - 2 n d - 1 (T (^§±)) 

r(¥)rm r /d-5\ r /7-dX 



£ ) " (7* • n - 6)^ .,±L 2 r( E^ r (±) r (§ 



7 3 6 4 647r3(r(|)) 2 r (f ) 

9aV r(|) 2 



F„J 320V2n 7 / 2 



Mlml (7fi • r 2 - 6) (<f a ■ cf 2 ) . (28) 
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III. RESULTS 



By repeating the analysis of the previous section, one can show that only the diagrams 
displayed in the Appendix produce extra contributions when non-relativistic propagators are 
used instead of static ones. The contribution of some of the diagrams containing two local 
vertices (from (jBJ)) can be extracted from [59J. In particular we have checked that the ratio 
of the pieces proportional to fractional powers of the masses arising from the diagrams Fig. 
1 b) and Fig. 1 f) in that reference agrees with ratio of (A. 13) and (A. 16) in the Appendix, 
as it should. There are further diagrams involving two-pion vertices with contributions in 
this region, but they are suppressed by powers of m^/M^ with respect to the ones displayed 
in the Appendix. This is due to the fact that two-pion vertices go with time derivatives 



(~ m^) rather than with space ones (~ i/m^Miv). Adding up all remaining contributions, 
and using the identities 

N^T a aNN^r a aN = -3N*NN*N 

N^T a NN^r a N = —2N^NN^N — N^erNN^crN (29) 

we obtain the following contribution to the (momentum space) potential 

V = 10^/2 M " m " { \{y) 6 + (§r) 2 \°s - C ^ ( 2 (* • **) + 4 ) + °t (2(* • 30 + 23 )] 



(30) 



This amounts to the following redefinition of Cs and Ct 

'3\2 . f-i /_ \ 6 / \2 



Cs - ^ Cs + w& M|ml { Kf ) 6 + (f ) 2 {cl - 4CsCt + 2Xi) 

= C s + ACs (31) 



Ct — Ct = C T + ^^ Mjml 1 2 (^j {-C S C T + C 2 T ) 

= C t + ACt (32) 

Hence the previous extractions from data of these coefficients must be corrected according 



to the formula above 



64|]. For instance, we may use the extraction found in ref. [7]. There 



we may find different results for Cs and Ct in the np channel coming from a fit done for 
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different values of the cut-offs A and A, which enter the Lippmann-Schwinger equation and 
the spectral-function representation of the two-pion exchange potential, respectively. The 
approximate results are shown in Table I. 



LEC 


{450, 500} 


{600, 600} 


{450, 700} 


{600, 700} 


C S 


-10.7 


8.9 


-12.1 


3.4 


Ct 


-1.2 


5.3 


-0.6 


2.5 



TABLE I: The approximate values of the S-wave LECs Cs and Ct in the np channel at N 3 LO 
for the different cut-off combinations {A[MeV], A[MeV]}, from the fit in [7]. The values of the 
constants are in 10 -5 MeV -2 . 



If we use these values of Cs and Ct to find an estimate of the size of the corrections ACs 
and ACt, we find the values displayed in Table II. Clearly these are no small corrections. 



A 


{450, 500} 


{600, 600} 


{450, 700} 


{600, 700} 


AC S 


21.2 


54.1 


23.3 


23.5 


AC T 


-1.6 


-2.8 


-1.0 


-0.3 



TABLE II: The values of the corrections ACs and ACt calculated using Cs and Ct from the 
former table as input. We use g A = 1-29, F n = 92.4 MeV, M N = 939 MeV, and m w = 139 MeV. 
The values of the constants are in 10 -5 MeV -2 . 



1, 



In fact, if using the data from [7j] we solve equations (|31~|) and (1321) exactly, we find complex 
values for Cs and Ct in all cases. This would indicate the need to redo the fits in [7] taking 
into account these new contributions. 



IV. DISCUSSION 

At first sight the new contributions to the potential we have found may look irrelevant 
since they amount to redefinitions of local counterterms. This is probably so as far as the 
description of scattering data is concerned. However, they may be of practical importance 
at least in the following two issues: (i) they shift the values of Cs and Ct extracted from 
data, which is important in order to check the consistency of a given counting scheme, and 
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(ii) they indicate that there will be contributions going like m^ 4 , m q being the light quark 
masses, which should be taken into account if one aims at a precision calculation of the 



nucleon-nucleon scattering lengths from the lattice QCD using chiral extrapolations 



60]. 



Notice that the scale ^Jitl^Mn coincides (up to corrections tji^/Mn) with the minimum 
momentum of the nucleons to produce a pion at rest. This coincidence suggests the possibility 
of an alternative derivation of our results, by computing the three-body spectral function 
and then making use of a dispersion relation. 

In higher loop calculations, subdiagrams with nucleon energies scaling as m n and nu- 



cleon three-momentum as ^/m^M^ will appear. It is not difficult to convince oneself that 
such subdiagrams are produced by adding an extra one pion exchange or a contact term 
to the two-loop diagrams we have calculated. These additions amount to a suppression 
of a factor M^ 2 mJ 2 / A5L which, parametrically, is equivalent to m l J 2 /M l J 2 if A x ~ M^- 
However, in practice takes a value ~ l$irF%/g\ in these contributions, which makes the 
ratio M 3 J 2 ml /2 /A 2 x ~ 1.2. Hence, one may consider a modified power counting in which 
M^ 2 ml/ 2 /A£ ~ 1, so that all these higher loop calculations should be consistently summed 
up. Whether this is feasible or not and whether this might provide an explanation to the 
unnatural size of the scattering lengths will be left for future work. 

Analogous contributions, proportional to fractional powers of the nucleon mass, are also 
expected in models with (arbitrary) boson exchanges at two loops (and beyond), as recently 



considered in 



6 II ]. In this reference non-relativistic propagators are used to calculate the 



contributions to the amplitude coming from the iterations of one-boson-exchange and two- 
boson-exchange potentials, but not in the calculation of the three-boson-exchange potential 
itself, for which the static approximation is used. We believe that, if non-relativistic rather 
than static propagators were used for the latter, contributions with the above mentioned 
fractional powers of the nucleon mass would also be obtained. 



V. CONCLUSIONS 



Using the method of the threshold expansions, we have found contributions to the nucleon- 
nucleon potential which are missed if static rather than non-relativistic nucleon propagators 
are used in the calculation and have calculated the leading contributions of them, which ap- 
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pear at two loops. They produce large contact terms with a peculiar non-analytic dependence 
on the light quark masses. 
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APPENDIX A: RESULTS FOR INDIVIDUAL DIAGRAMS 
1. Diagrams with no contact terms 



— i 



9A 



r(!) s 



320 v / 2tt 7 /2 



M^mliflrx ■ T 2 - 6)(<7i • <7 2 ) 



(Al) 




r(!) : 



9a 



F^J 640^7/2 



M^ml{fx ■ f 2 + 6)(<7i • <r 2 ) 
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Diagrams with one contact term 



T 



= i 



r(I) 1 



9a 



FJ 240^^/2 



Mlml{2r 1 ■ f 2 - 3) 



x (C s (2*i • <r 2 - 3) - C T (7ai • a 2 - 6)) 




—i 



9a 



r(!) : 



i 7 ;/ 240V2 7T 7 / 2 



M£m|(2fi • r 2 + 3) 



x ((75(2,7! • <r 2 + 3) - C T (^i • <? 2 - 6)) 
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9A 



4 ■ < " 



r(!)' 



/ | — i - - _ _._ M£m% 



x [(3f! • r 2 ) (3C 5 (ai • B 2 ) + C T {2B 1 ■ a 2 - 3)) 

+(2fi • r 2 - 3) (C 5 (2a! • <? 2 - 3) + • a 2 - 6))] (A9) 



1 480^^/2 ' 



x [(fi • f 2 ) (C s (ai • <7 2 ) - C r (2ai • a 2 + 9)) 

+(2fi • r 2 + 3) ((75(2,7! • <7 2 + 3) - C r (*i • a 2 + 6))] (A10) 



1 I i|) 160^/2 M ^ m ' (fl ' ' ^) " C T (2a 1 ■ a 2 - 3)) 

(AH) 



(A12) 
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3. Diagrams with two contact terms 




—i 



9A 



Fj 20V2tt 7 /2 



x (C|5i • <r 2 - 2C S C T (28 X ■ & 2 ~ 3) + C£(7ai • a 2 - 6)) 



(A13) 




= 2 



r(!) ; 



i 7 ;; 40V2V/2 



x [3(3C| + 2C 5 C T (<? 1 • «t 2 ) + (71(2,7! • <r 2 - 3)) 
+(ri • r 2 )(C 2 s (^i ■ <r 2 ) + 6C 5 C T - CUffi ■ °2 - 6))] 




9A 



r( 



!) 2 



i 7 ;; 20v / 2tt 7 /2' 



(A14) 



x [3(3C| - 2C s C T (<7i • <t 2 ) + C|(2ai • <t 2 + 9)) 

+(ri • f 2 ) (CK?! • <r 2 ) + 2C S C T {23 1 ■ & 2 + 3) - CK^! • a 2 + 6))] 

(A15) 




—2 



9r(f) 2 

40^^/2 



M*m*(Cj + 2<7 s C r (a 1 • * 2 ) - C^(2a 1( x 2 - 3)) 



(A16) 
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[63] An analogous situation has been identified in the contexts of hadronic atoms 43J and heavy 
quarkonium [441 ] . 

[64] There are further contributions from the iteration of the potentials and from higher order 
potentials themselves which eventually amount to redefinitions of Cs and Ct- These contribu- 
tions, however, have already been taken into account (see, for instance, [3]) in the fits to data. 
We restrict our discussion to the new type of contributions found in this work. 
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